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Abstract 

The phase separation process which follows a sudden quench inside the co- 
existence region is considered for a binary fluid subjected to an applied shear 
flow. This issue is studied in the framework of the convection-diffusion equa- 
tion based on a Ginzburg-Landau free energy functional in the approximation 
scheme introduced by Bray and Humayun [Phys. Rev. Lett. 68, 1559, (1992)]. 
After an early stage where domains form and shear effects become effective 
the system enters a scaling regime where the typical domains sizes Ln, Lj_ 
along the flow and perpendicular to it grow as i 5 / 4 and i 1//4 . The structure 
factor is characterized by the existence of four peaks, similarly to previous 
theoretical and experimental observations, and by exponential tails at large 
wavevectors. 
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I. INTRODUCTION 



The process of phase-separation occurring in binary systems quenched below the co- 
existence line has been matter of intensive investigations since many years. The overall 
phenomenology is nowadays reasonably well understood. In an early stage, whose prop- 
erties depend on the system specific details, domains of the two phases form; the relative 
concentration of the two species well inside the domains quickly saturates to a value very 
similar to the equilibrium composition. Then a crossover leads to the late stage dynam- 
ics which is characterized by the growth of ordered regions of typical size L(t) and by a 
statistically invariant morphology. This phenomenology is at the basis of the well known 
scaling property according to which a single length L(t) is asymptotically relevant: system 
configurations at different times are simply related by a scale transformation with a scale 
factor L(t). In particular, for the equal time density-density correlation function one has 

(1) 

meaning precisely that G is invariant as a function of scaled lengths r/L(t). The growth 
of L(t) is of the power law type L(t) ~ t 1 ^. The dynamic exponent z depends on the 
specific microscopic mechanisms responsible for the segregation process; in general different 
dynamical regimes can be observed, characterized by different values of z, by changing the 
experimental time window fl]]. When the system enters the late stage the memory of the 
initial state is lost and, due to the presence of one single growing length, its properties become 
universal. Then the dynamical exponents, the scaling function g and other observables 
become independent on the system specific details and are only determined by a few relevant 
parameters. Besides the spatial dimensionality, the presence of conservation laws is a relevant 
feature. For binary fluids, to which we restrict our attention, the concentration of the species 
is conserved as the phase-separation proceeds. The initial stage of the scaling regime, before 
hydrodynamic effects become relevant, is the subject of this paper. For viscous fluid this time 
window can be sufficiently wide to be detected experimentally. In this condition the phase- 



G(r,t)=g 



L(t) 
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separation process is ruled by the evaporation, diffusion and subsequent recondensation of 
monomers |2j and the exponent z is known to be 3, independent on the space dimension d. 

The spatial properties are encoded into the scaling function g(x) or, equivalently, into 
its Fourier transform f(q). For large q, f has a power-law tail of the form This is 

the famous Porod's law |3[], long recognized as arising from sharp domains walls. 

Another relevant parameter affecting the scaling properties is the number of components 
N of the order parameter 0. For binary fluids is a scalar field. Most of the analytical 
techniques developed to study the kinetics, among which the Bray-Humayun scheme that is 
considered throughout this article, are however suited for vectorial systems with large- N. Let 
us briefly discuss the main features of vectorial systems that will be recalled in this paper. 
For systems with a continuous symmetry, rotations of the order parameter on the O(N) 
symmetric manifold of the local ground state of the Hamiltonian, instead of evaporation 
and condensation, is the important growth mechanism. This changes z to 4, regardless 
both of d and N (apart from exceptional cases Concerning the large-g behavior of / 
a generalization of the Porod's law, f(q) ~ q~( d + N ) is obeyed for iV < d. This form 
relies on the existence of stable localized topological defects. Much less is known about the 
tails of / when these defects are not present, for N > d. Numerical simulations found a 
stretched-exponential decay f(q) ~ exp[— q v \. 

Despite this rather simple phenomenology the problem of deriving the scaling properties 
from first principles is still open. On the analytical side, the only soluble model || for 
conserved order parameter is the one with N = oo. When considering this unphysical limit 
one hopes that the gross features of the kinetic process are retained for N = oo leaving 
further refinements to 1/N perturbation theories. However, although the global behavior of 
the N = oo model is qualitatively resemblant [7[ of real systems, the scaling symmetry ([[]) 
is violated already at a qualitative level. The explicit solution of the model || shows in fact 
that instead of the standard scaling ([]]) one has a multiscaling symmetry produced by the 
existence of two logarithmically different lengths. Actually the ordering process with N = oo 
is more a condensation in momentum space ||, than a genuine phase-separation process. 
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These results shadow the possibility to develop a successful theory of the segregation process 
by systematic 1/N expansions around the large- N limit. 

In this scenario a prominent role is played by a class of approximate theories known as 
gaussian auxiliary field theories, originally introduced by Mazenko |J . The essence of these 
is a non-linear mapping between the order parameter <j)(x,t) and an auxiliary field m(x,t). 
The actual form of the relationship m(4>) is largely arbitrary and, depending on different 
physically motivated choices, different schemes have been developed jlOfl . The basic idea is 



that with a proper choice of m(4>) the equation obeyed by the auxiliary field can be treated 
to lowest order in a mean-field approximation while the basic non-linearities of the evolution 
are still retained through the non-linear mapping m(0). This ad hoc assumption is generally 
uncontrolled and its use is justified a posteriori by the satisfactory quality of the predictions. 
Starting from the method of Mazenko, Bray and Humayun (BH) have derived ||TT] 



a non- 



linear closed equation of motion for the correlation function G(r, t) within the framework 
of the 1/N perturbative expansion. Letting iV — > oo in the BH equation one recovers 
the large-iV model previously discussed; the BH approach, however, does not amount to 
a straightforward 1/iV-expansion because all orders in 1/N are reintroduced through the 
gaussian auxiliary field approximation. On one side, the uncontrolled character is the limit 
of the BH scheme, since it is not clear how to improve the approximation, but on the other 
side, represents a powerful tool to overcome the shortcomings of standard perturbation 
theory previously discussed. As shown by BH their equation reinstates standard scaling ([[]) 
as a truly asymptotic symmetry, leaving to multiscaling a merely pre-asymptotic role. An 
additional feature at variance with the N = oo limit is the form of the scaling function f(q), 
which, in the BH approximation falls exponentially at large momenta, as opposed to the 
quartic exponential decay of the N = oo model. 

In this paper we are interested in the diffusive phase separation process of a binary fluid 
in the presence of an applied shear flow, namely an imposed fluid velocity along x with a 
constant gradient along y. The presence of the flow radically changes the behavior and even 
the phenomenology is much less understood. A first trivial effect is the anisotropic defor- 



mation of the growing pattern that appears greatly elongated along x JT2|,|T3|,[Tj,|T5[ . Then, 
when dealing with growth laws one has to specify the length of the equilibrated regions 
for each direction, L x (t), L y (t), L z (t), and in general L x (t) » L y (t), L z (t). A more subtle 
point, that may have important consequences, is the following: stretching of the domains 



causes ruptures of the network [H that may render the segregation incomplete. Actually, 
the characteristic lengths in the different directions could keep on growing indefinitely, as 
without shear, or due to domains break-up the system may eventually enter a stationary 
state characterized by domains with a finite thickness. In general there is no agreement on 
this point both from the experimental Jl5j and theoretical point of view. On the theoreti- 



cal side large- N calculations [fT6|JT7[] show the existence of a multiscaling regime with ever 



growing lengths in all directions. However numerical simulations |18| cannot still establish a 
clear evidence due to finite size and discretization effects. Another intriguing feature which 
is probably related to the stretching and break-up of the domains is the observation of an 
oscillatory pattern in the observables. Let us consider the excess viscosity At] as an exam- 
ple. Stretching of domains require work against surface tension leading to an increase of 
Ar). However when a domain breaks up the stress is released since the fragments are more 
isotropic; this produces a decrease of At]. If repeated in time this mechanism produces an 
oscillatory pattern that is observed in simulations [18| and possibly in experiments fll9 |. In 



the large- N limit it can be shown that this process has a periodicity in logt. Similar pat- 
terns are also reported in apparently different contexts as fractures in heterogeneous solids 
| 2Cfl or in stock market indices ||21|| , suggesting a possible common underlying interpretation. 



However we do not have nowadays analytical tools to enlighten this feature. 

The spatial properties of the mixture are encoded into the correlation function G(f, t) or, 
equivalently, into its Fourier transform C(k,t), the structure factor. The latter is directly 
measured with scattering techniques and is more suited for comparison with the experiments. 
The form and the evolution of C are themselves modified by the presence of the flow. When 
the shear becomes effective, after an initial isotropic regime where C has the shape of a 
circular volcano, the tip of C is deformed into an ellipses. This is expected as a consequence 



of the anisotropy, and is generally observed in experiments. In some cases the presence 



of four pronounced peaks on the edge of C can be detected p2fl. The existence of these 



multiple maxima has been also observed in numerical studies [18] of phase separation and 



in different physical systems such as stationary microemulsion in shear flow [2^] , suggesting 
that this is a rather general feature of sheared fluids. A peak in the structure factor is 
generally interpreted as the signature of a characteristic length proportional to the inverse 
of the wavevector where the peak is located. In this anisotropic case to each peak one 
associates one typical length for each space direction. Taking into account the symmetry 
k — > —k this suggests that a couple of characteristic lengths for each direction is present: It 
is not completely clear, however, to which extent this interpretation can be pushed. 

In this partly unclear scenario the development of theoretical tools for reliable predictions 
to be tested in experiments is important. In the case without shear the BH scheme turns 
out to be a reference theory capable to describe the main ingredients of the phenomenology. 
Therefore in this paper we address a numerical solution of the BH equation in the presence 
of shear flow. We will compare its behavior with previous results pertaining to the large- N 
limit or to numerical simulations of the fully non-linear scalar model. Our results confirm 
the occurrence of a scaling regime with growing lengths in every direction. The dynamical 
exponents are found to coincide with those of the N = oo model. The four peak pattern 
for the structure factor is exhibited, similarly to what observed both for N = oo and in 
the scalar model, the actual form of C being, however, rather different in the three cases. 
Concerning the oscillating behavior of the typical lengths and of the excess viscosity, on the 
other hand, the BH model behaves quite differently, in that these are strongly inhibited or 
even absent, at variance with N = oo and N = 1. 

This paper is divided in 4 Sections: In Sec. [TT] we introduce the BH equation and discuss 
some numerical details about its integration; in Sec. [TT1] the results of the numerical solution 
are presented and a comparison with different approaches is proposed. In Sec. |TV] we draw 
the conclusions of this work. 
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II. THE MODEL 



We consider a system described by an order parameter with iV components. The evolu- 
tion is described by the convection-diffusion equation [JT| that generalizes the Cahn-Hilliard- 
Cook equation in the case of an applied velocity field v 

« + v ■ w , ms, *)] = rv^qM + , ( f , t) (2) 

where <fr = {4>a} is the vector order parameter that in the scalar case represents the concen- 
tration difference between the two species of the fluid, x = {x, y, z] is the space coordinate, 
T is a transport coefficient and r\ is a gaussian stochastic field with expectations 

(r ] {x,t)) = 

(t](x, t)rj(x, t')) = -2TTV 2 5(x - x)8(t - t') (3) 

describing thermal fluctuations. The symbol (. . .) means an ensemble average. T is the 
equilibrium Ginzburg-Landau free energy 

nm = J dx + jl + (4) 

where 2 = J2 a ( t ) a ( t ) a- The parameter r distinguishes between the mixed state with r > r c (T) 
(r c (T) < 0) and the phase separated states with r < r c (T). For a plane shear flow the velocity 
term is given by 

v = jye x (5) 

where e x is the unit vector in the flow direction x and 7 is the shear rate. 

In the approximation framework developed by BH 0] the following equation of motion 
can be derived for the equal time correlation function G(r, t) = {4> a {x, t)4> a (x + r, t)) 

1 



dG(f, t) + ^ dG(f,t) = _ 2y2 , _ 



G(f,t) + -G 3 (f,t) 



(6) 



dt dx 

where we have dropped the component indices due to internal symmetry. In deriving Eq. (^) 
from Eq. (||]) we have let T = 1 and u — 1 since this simply amounts to a redefinition of 



time and space scales and we have neglected the thermal disturbance i] because it can be 



shown that the temperature is an irrelevant parameter below the coexistence line. 

The quantity R(t) is a function of time that is asymptotically determined by the requirement 
lim^oo G(0, t) = S eq = —r/g. Here we use the self-consistent determination of R(t) 

R(t) = g[G(0,t)-S eq ] (7) 

that is usually adopted [|TT1 , ^5| . For N = oo, namely dropping the cubic term on the r.h.s., 



Eq. (^[) becomes the N = oo model. 

We have considered the evolution of Eq. (||) with N — 1. Some comments about the role 
of N are in order. Eq. (|6]) is derived in an 1/N perturbative framework. Then, in principle, 
the quality of the approximation is expected to be better for large N and, in any case, this 
whole approach is meaningful for vectorial systems with continuous symmetry. Given that 
the nature of the symmetry is relevant for the scaling properties, as discussed in Sec. |I[ one 
cannot pretend to push this scheme down to the physically relevant case with iV = 1 in a 
completely successful way. However, besides numerical simulation that are quite difficult 
and not particularly instructive in this case, most of the understanding of phase separation 
with shear comes nowadays from the N = oo model, which surely suffers of much more 
profound shortcomings than the present approach, if used to infer the properties of physical 
systems. The BH scheme represents a first step beyond N = oo, despite of course being 
not conclusive. It must be recalled, moreover, that at least without shear the BH equation 
substantially reproduces the behavior of the system with N = oo, including multiscaling, in 
a preasymptotic time domain t < t*(N) ~ (AiV) 4 / d (hi N) 3 after which the cubic correction 



in Eq. (jty becomes effective pT| . The choice N = 1, besides being appropriate for physical 
systems, is also the most efficient in order to amplify as much as possible the role of finite- N 
corrections, given the limited timescale numerically accessible. 

Eq. (^|) have been solved on a two dimensional lattice of size L = 4096, the lattice 
constant being Ax = 1, via a finite difference first order Euler scheme with At = 0.01. 
Periodic boundary conditions have been adopted. The initial condition G{r, 0) = A<5(f), A 
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being a constant, corresponds to a quench from an equilibrium state at infinite temperature. 
The parameter of the free energy we use are — r = g = 1; different choices correspond to 
a redefinition of the order parameter amplitude. We will present data relative to the case 
7 = 0.01 and A = 0.01; we have tested that other choices produce similar results. 
The structure factor C(k,t) = (<j>(k,t)(f>(—k,t)) obeys the equation 

where D(k,t) is the Fourier transform of G 3 (f,t). Due to the structure of this Equation, 
however, it is more efficient to compute G from Eq. ([]) and to obtain C by Fourier transform. 
From the knowledge of C one can extract the characteristic lengths along x and y as 

' J dk\k x \C{k,t) 



L x {t) 



JdkC{k,t) 



JdkC(k,t) 

and rheological indicators such as the shear stress 



L y {t) = I / dfc IWM) , (g) 



dk 

°*y(t)= I J^yk x k y C(k,t) (10) 



For a steady flow the excess viscosity is defined in terms of the shear stress as A^(t) 
--f~ l a xy (t). 



III. RESULTS 

The anisotropic character of the evolution, with the ordered regions aligning along the 
flow direction, is reflected in the behavior of the correlation function G, shown in Fig. ([[]). 
G has a peak at r = and is continuously stretched along the x-direction, assumining 
increasingly eccentric elliptical patterns. 

Snapshots of the structure factor at the same times are shown in the k x , k y plane in 
Fig. (0). Initially, not shown here, C develops the form of a circular volcano, as in the 
case without shear. A later observation, at jt — 1, shows a deep in the volcano profile 
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along the direction k y = —k x . This deep becomes progressively more pronounced until C 
appears separated into two symmetric foils. At times of order jt = 5 in each foil a couple of 
maxima are built up. We denote by the letter A the peak that in the figure at time jt = 5 
is higher. With respect to the other maximum, denoted by B, A is characterized by having 
\k x \ » \k y \, whereas in B \k x \ ~ \k y \. The relative height of these maxima changes in time. 
This is observable in Fig. (0) at time jt = 9: now B is prevailing. 

Up to this time the qualitative evolution of the structure factor is comparable to what 
reported for N = oo |17| and for the numerical solution [18| of the full model equation with 



N — 1. Computer limitations prevent the observation of the system on longer times. Only 
in the N = oo case the dynamics can be followed by numerical integration of the equations 
on larger timescales and it can be shown that the repeated prevalence of either A or B 
maxima continues cyclically in time. It turns also out that the recurrent dominance of A 
and B is periodic in the logarithm of the strain r yt. A physical interpretation of this pattern 



in terms of stretching and break-up of domains is proposed in |L8 
The spherical average 

C sf (k, t) = AT 1 J d6C{k, t), (11) 

where tan# = k y /k x and M is the number of lattice points contained in a circular shell of 
width 7r/256 centered around k, and the averages along k x or k y 

C x (k y ,t) = L- 1 J dk x C(k,t) (12) 

and 

C y {k x , t) = L- 1 J dk y C{k, t) (13) 
are plotted in Fig. (Rf). This figure shows that in the large-fc tail the structure factor decays 



exponentially, as already observed without imposed flow ||25|| , at variance with the faster de- 
cay observed for N = oo [|T^]. As already mentioned, the absence of power-law tails is related 
to the absence of stable localized topological defects in the Bray-Humayun approximation. 
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Next we consider the growth laws. For N = oo, in the large time region fl6|J17| L x and 
L y keep increasing with (logarithmically corrected) power laws L x ~ 7(t 5 / hit) 1 / 4 , L y ~ 
(t/lnt) 1//4 modulated by a log-periodic oscillation. The short time part of this pattern is 
observed in Fig. (f|b). The same growth exponents are expected for the BH model (but 
without logarithmic corrections); this is because the model amounts to an 1/N perturbative 
expansion and, therefore, it is accurate for vectorial systems whose growth exponents belong 
to the same universality class of N = oo. Actually the value of the exponents can be easily 
determined by means of a scaling analysis. If standard scaling holds the structure factor can 
be cast as 



C(k,t)=L x (t)L y (t)f(X,Y) 



(14) 



where X = k x L x (t) and Y = k y L y (t) and / is a scaling function. From ([TJ]) one also has 



D(k,t) =bL x (t)L y (t)h(X,Y) 



(15) 



where b is a constant and h another scaling function. Inserting expressions (114 , 15]) into 
the equation of motion of C (H) and assuming that asymptotically L x prevails over L y one 
obtains 



2Y A f(X,Y) = -2R(t)L y (tf 



1 dL y (t) 4 
4 dt 



1 + 



f(X,Y) + -h(X,Y) 

d[lnL x (t)] 
dt 



V 2 | L y (tf df(X,Y) 



d[\nL y {t)\ 
dt 



(16) 



The l.h.s. of Eq. ([16]) does not depend explicitly either on time and 7. Then in order to 
fulfill this equation for each X, Y both t and 7 must drop out from the r.h.s. Then one has 

L x {t) ~ 7*i 

Ly(t) ~ t\ (17) 

R(t) ~ rl 

As already anticipated the exponents are those pertaining to systems with continuous sym- 
metry. Furthermore one can observe that the shear rate enters linearly in L x but the growth 
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in the shear direction y is unaffected by 7. For completeness we recall that in the full 
scalar model one expects different exponents, namely L x ~ 7^ 4//3 , L y ~ t 1 / 3 from scaling 
p~6|] or renormalization group |18| arguments, but these are not yet clearly observed in sim- 
ulations. Moreover the very existence of an asymptotic scaling dynamics is debated in this 
case. Fig. (^a) shows that initially both L x and L y start growing with an exponent which is 
consistent with the value 1/4 as in absence of shear. In this time regime the fluid essentially 
does not feel the presence of the flow, as already discussed for the structure factor. Then, 
starting from ~ 0.5, L y keeps growing unexpectedly larger than L x untill at yt ~ 2 it 
decreases. Up to this point the BH model roughly resembles the N = 00 case, despite the 
fact that the overshoot of L y is much more pronounced in the latter case. Fig. (|4]c), instead, 
shows that the full scalar model behaves quite differently, with a pronounced decrease of L y 
in the range of times where in the other models a fast growth was observed. For longer times 
the asymptotic stage is entered. The numerical integration of the BH model presented inso- 
far shows that Eq. fli7| ) is obeyed but, at variance with the other two cases, the oscillatory 
pattern is absent or, at least, strongly depressed. 

Now we turn to the discussion of Fig. (^) where the behavior of the excess viscosity is 
reported. In all the cases considered Arj initially grows to a maximum, then falls negative and 
raises to a second maximum. This pattern has been observed also in experiments [|IIJ and is 
referred to as double overshoot. After the second maximum the asymptotic stage is entered. 
On the basis of scaling arguments one expects that A77 scales as 7 -1 times the inverse of 
the typical volume L x L y of the domains of the two species. Then Ar/ ~ 7~ 2 £~ 3 / 2 for vector 
order parameter and Ar/ ~ 7~ 2 £~ 5 / 3 for the scalar model. The former exponent is observed in 
Fig. ([|a) for the BH model, although in the very last part of the simulation. For N = 00 the 
same exponent (with logarithmic corrections) has been exactly computed |16[ and observed 
numerically on larger timescales, but modulated by log-time periodic oscillations |I7fl . In 
the scalar case a definite exponent cannot be extracted from simulations, but the oscillations 
are observed. 



12 



IV. CONCLUSIONS 



In this Article we have considered the behavior of a binary fluid quenched below the 
mixing temperature in the presence of an applied shear flow. The description of this system 
has been carried out in the context of the continuum convection-diffusion equation based on 
the Ginzburg-Landau equilibrium free energy. The nature of our approach, which neglects 
hydrodynamic effects, limits the domain of applicability of the model to the diffusive regime 
which takes place at short times. We have studied this model in the approximation scheme 
introduced by Bray and Humayun. This amounts to an expansion in 1/JV on top of the 
auxiliary field theory approximation originally introduced by Mazenko. A comparison is 
also presented with two different approaches: the N = oo limit and the full scalar model, 
which is the most appropriate for fluids. The first has the great advantage to be fully sol- 
uble analytically, al least in the large time domain, providing an unambiguous description. 
Furthermore it can be easily studied numerically from the instant of the quench onwards. 
However, besides the fact that the case N = oo is far from physical, it has also proved 
to violate the qualitative symmetry of standard scaling in favor of multiscaling, rendering 
this model insecure as an approximation to real systems. On the other hand the full scalar 
model is well suited for describing fluids but it cannot be studied analytically. Numerical 
simulations, although important, have not yet provided clear evidences of the basic phe- 
nomenology, due to technical limitations. The BH model is ideally located in between these 
two: It shares with the large- N limit the property of being a closed equation for an ob- 
servable, the correlation function, that is already an ensemble average. Then, even if one 
resorts in the end to a numerical solution, one still has the great advantage to avoid the 
average over a large number of realizations. Secondly, due to its perturbative character, it 
is believed to represent a step forward with respect to the case N = oo. This statement has 
been proven to be correct at least without shear. 

Given the essence of the BH scheme, one does not expect to obtain a fully reliable 
approximation to real systems with N — 1. The spirit of this work is more to establish 
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if and to which extent its phenomenology compares better to the numerical simulations 
of the physical model than to the case N = oo. From the analysis of the behavior of the 
characteristic lengths and of the excess viscosity we can conclude that the BH model behaves 
more like the N = oo case with damped oscillations than like the full scalar model, despite 
that, with the choice N = 1, the corrections to the case N = oo contained in Eq. are 
effective almost from the beginning, as it is clear from Figs. (@,^). This is true not only 
for the value of the dynamic exponents, but also for the qualitative behavior of the typical 
lengths. 

On the other hand, regarding the oscillating character of most observables, the BH 
model behaves rather differently from both N = oo and iV = 1 in that strong oscillations 
are observed in the latter two while they are practically absent in the former. Under this 
respect the BH model is atypical and cannot be regarded as a bridge between N = oo and 
N = 1. It is possible that the BH scheme resembles more the behavior of vectorial system 
with a finite N under shear, but this is a totally unexplored field. 

Regarding the important issue of the existence of an asymptotic scaling regime charac- 
terized by ever growing lengths L x ,L y , in the time range accessed in this work we do not 
see any saturation effect and both L x and L y keep firmly growing with the expected power 
laws. In the BH scheme, then, the presence of an asymptotic scaling regime is confirmed 
while a stationary state with domains of a finite thickness is not observed. 
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FIGURES 
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FIG. 1. The real-space correlation function G(r,t) is shown (only a part of the whole lattice) 
at four different values of the strain, jt = 1, 5, 7, 9. x is on the horizontal axis. 
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FIG. 2. The structure factor C(k, t) is shown at four different values of the strain, jt = 1, 5, 7, 9. 
k x , on the horizontal axis, increases from left to right; ky increases toward the back of the figure. 
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FIG. 3. The structure factors averaged along the k x and k y directions and spherically are shown 
for yt = 10 (•), 15 (*), 20 (A), 25 (o) and 30 (*). 
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FIG. 4. The evolution of the characteristic lengths L x (t), L y (t) is shown for a) the BH model, 
b) the N = oo limit, c) the full scalar model (numerical simulation). In a) and b) the straight lines 
represent the expected power law behaviors L x (t) ~ i 5 / 4 , L y (t) ~ t 1//4 . 
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FIG. 5. The evolution of the excess viscosity Arj(t) is shown for a) the BH model, b) the N = oo 
limit, c) the full scalar model (numerical simulation). The straight lines represent the power law 
behavior Arj(t) ~ i~ 3 / 2 . 
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